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Abstract 

Systematic properties of the first excited state are presented for various ring 
sizes and spin quantum numbers which are only partly covered by the theo- 
rem of Lieb, Schultz and Mattis. For odd ring sizes the first excited energy 
eigenvalue shows unexpected degeneracy and related shift quantum numbers. 
As a byproduct the ground state energy as well as the energy of the first 
excited state of infinite chains are calculated by extrapolating the properties 
of only a few, finite, antiferromagnetically coupled Heisenberg rings using the 
powerful Levin sequence acceleration method. 
PACS: 75.10.Jm, 75.40.Cx 

I. INTRODUCTION AND SUMMARY 

Exact diagonalization methods make it possible to investigate small spin rings for 
various numbers N of spin sites and spin quantum numbers s, for instance in the Heisenberg 
model. The symmetries of the isotropic Heisenberg Hamilton operator allow to decompose 
the Hilbert space H into a set of mutually orthogonal subspaces H.(S, M, k) according to 
the quantum numbers of the total spin S, the total magnetic quantum number M and the 
translational quantum number k of the cyclic shift operator. Since these subspaces are much 
smaller, a complete or partial diagonalization of the Hamilton matrix is feasible. 

The interest in the Heisenberg model, which is known already for a long time, was 
renewed by the successful synthesis of small magnetic molecules, among them nearly perfect 
ring structures of paramagnetic ions like Fe 3+ One can say that in the majority of these 

molecules the ions couple antiferromagnetically and the spectrum is rather well described 
by the Heisenberg model with isotropic next neighbour interaction. 

Looking at the properties of the calculated spin rings one realizes that not only the ground 
states, but also the first excited states share systematic properties which are understood 
only for Heisenberg spin rings of an even number of spin sites (bipartite systems), i.e. they 
can in part be derived from the theorem of Lieb, Schultz and Mattis PJlOfl. Non-bipartite 
rings which have an odd number of sites, and thus can be characterized as frustrated, show 
unexpected properties in degeneracy and translational quantum number. 

Knowing systematic rules for quantum numbers also of non-bipartite systems would be 
very useful for comparison of theoretical results with measurements in ESR/EPR, torque 
magnetometry or neutron scattering, see e.g. 0,01 • One could employ knowledge about 
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quantum numbers of ground and first excited states in order to understand the thermal 
behaviour of quantities like the magnetic susceptibility. In addition these exact values may 
help to improve low temperature approximations. Usually the high temperature behaviour 
of observables is well known, e.g. from classical spin dynamics |L2|, but at low temperature 
such approximations are poor. The knowledge of ground and first excited states could 
already be sufficient for a considerable improvement. 

Having evaluated the spectra of small Heisenberg rings with isotropic next-neighbour 
interaction one can approximate the infinite chain limit, which for the s — | ground state 
is known as the Bethe-Hulthen limit [13,13]. Because the sequences converge rather slowly 
and only a limited number of energy eigenvalues can be evaluated, the Levin u-sequence 
acceleration method |I5| , |T6| is employed, which leads to impressive estimates of the antifer- 
romagnetic ground state energies as well as of the excitation gap for infinite rings or chains 
of larger spin quantum numbers. 



II. SYSTEMATIC PROPERTIES OF THE SPECTRUM 

The Hamilton operator of the Heisenberg model with antiferromagnetic, isotropic next 
neighbor interaction between spins of equal spin quantum number s is given by 

JV 

H = -2 J 10) + , Vx:s(x) = s, J<0, N + 1 = 1 . (1) 

2 = 1 

H commutes with the total spin S and its three-component 5* 3 . In addition H is invariant 
under cyclic shifts generated by the cyclic shift operator T. T is defined by its action on the 
product basis (|3]) 

T | mi, ... , mjv-i, m N ) = \ m N , mi,..., mjv-i ) , (2) 

where the product basis is constructed from single-particle eigenstates of all s 3 (x) 

s 3 (x) | mi, ... , m x , m N ) = m x \ m 1 , m x , m N ) . (3) 

The translational quantum number k = 0, . . . , N — 1 modulo iV labels the eigenvalues of T 
which are the iV-th roots of unity 

z = expj-z^j . (4) 

Exact diagonalization methods HQ allow to evaluate eigenvalues and eigenvectors of H for 
small spin rings of various numbers N of spin sites and spin quantum numbers s. Systematic 



investigations |17H2T| revealed interesting properties of ground state quantum numbers, 
compare table |, which only for bipartite rings can be explained by the theorem of Lieb, 
Schultz and Mattis [§,|Kj]. The ground state spin quantum number and the degeneracy, for 
example, depend solely on iV-s. If N-s is integer, then the ground state has S = and is 
non-degenerate; if N-s is half integer, then the ground state has 5=1/2 and is fourfold 



degenerate pi . 
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It appears that for the properties of the first excited state such rules do not hold in 
general, but only for "high enough" N (N > 5). Then, as can be anticipated from table |, 
we can conjecture that 

• if N is even, then the first excited state has S = 1 and is threefold degenerate, and 

• if N is odd and the single particle spin is half-integer, then the first excited state has 
S = 3/2 and is eightfold degenerate, whereas 

• if N is odd and the single particle spin is integer, then the first excited state has S = 1 
and is sixfold degenerate. 

Except for small odd N also the cyclic shift quantum numbers k of the first excited state 
show an interesting regularity. 

• For odd N > 7, k assumes a certain value for all integer spins and another value for all 
half-integer spins. We conjecture that the k quantum numbers for half-integer spins 
are k = 3 [(N + 1) /4J and k = N - 3 [(N + 1) /4J . [(N + 1) /4J symbolizes the greatest 
integer less or equal to (N + l)/4. For integer spins numerical data are poor but it 
seems that k is as close as possible to N/2, i.e. k = \_N/2\ and k = N — \_N/2\. 

• For even N the shift quantum number k is N/2, if N/2 is also even; if N/2 is odd, 
k = for half- integer spin and k = N/2 for integer spin. 



For spin-i-rings these properties may be also derived using the Bethe ansatz |iTi].[i~f 



III. INFINITE CHAIN LIMITS 



Besides the importance of the above presented results for magnetic molecules 
|I|j2^j7| j23|JTl|j2^1 , the obtained energy eigenvalues enable us to estimate the antiferromag- 



netic ground state energy E (N) in the large N limit for a variety of spin quantum numbers. 
Of course this calculation cannot compete with nowadays DMRG results, but reaches as- 
tonishingly close. 

As one can see in table | or Fig. || (l.h.s.) the convergence of energy eigenvalues with 
N is rather slow. Therefore, an improved estimate is calculated using the Levin w-sequence 



acceleration method [15], which is appropriate if the series elements form an approximately 
linear function in -4 with a certain positive k, which to first approximation is the case for 
the sequences of table III This observation is in accordance with the Wess-Zumino-Witten 



model, see e.g. [p5fl , which yields 

E ot 



Let us denote the elements of our series by U\, C/ 2 , U 3 , In order to construct monotonic 

sequences the alternating series for a certain spin are divided into two monotonic subseries 
e.g. for s = \: U x = E (N = 2), U 2 = E (N = 4), U 3 = E (N = 6), . . . . The differences 
between successive sequence elements are labeled U\ — Ui, Ui — Ui — U\, U3 — U3 — U2, 
Then the Levin u-estimate using n series elements reads 
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U\ n ] = ^ k=l[ ' \hll !* (6) 

The following example for the sequence built of the energies for even numbers of spin s = 1/2 
demonstrates how fast the Levin ^-estimate approaches the correct limit 

{U[n}} = {1.5, 1.2, 0.8735510038, 0.8885066176, 0.8858640679, 0.8863719562, (7) 
0.8862817068,0.8862961998,0.8862941347,. . .} -> 0.8862943611 = 21n2 - - . 

Roughly one can say that, if the deviations from a power-law behaviour are small, n sequence 
elements lead to an accuracy of about n — 4 digits, see Fig. |T] (r.h.s.). Of course the limit 
can only be as accurate as the individual sequence elements. 

The procedure is also applied to spin values s = 1, 3/2, 2, 5/2 and the results are depicted 
in Fig. 0. The hatched areas indicate the interval according to the gained accuracy of the 
sequence acceleration method, the attached number denotes the mean value of the interval. 
Since the sequences built from values for even N converge much faster, they determine the 
results. 

In the upper left of Fig. |2] the estimate for the infinite chain limit for spin s — 1 is shown. 
The estimated ground state energy of E /(NJ) = 2.802967 ± 0.0000005 agrees very well 
with the result 2.802968077942(4) found in Ref. Q and also with other DMRG and TMRG 
calculations, see e.g. |3T| , p2| . Also for s = 3/2 the limit is rather well approximated, the value 
of E /(NJ) = 5.65681 ± 0.00001 suggests new discussion of the results 5.65666 ± 0.00002 of 



Ref. and 5.658 of Ref. p2[. For these sequences previously obtained exact diagonalization 



results have been used, too [pB} - p9| , |3T 



The excitation energies of the first excited state, see table |, enable us to approach the 
gap for infinite chains or rings. It vanishes for half integer spins and remains finite for integer 
spins, Haldane conjecture P3LE5J. The following example shows the Levin u-sequence for 
s — 1/2 and even N 

{U[n}} = {4., 3., -2.464274955, 0.26749212, 0.05231106824, -0.04294415611, (8) 
0.01234278872, -0.003416797715, 0.0001536344576363, . . . } -»• . 

In Fig. | one can see that the behaviour of the sequence is much smoother for even N, 
whereas the somewhat strange behaviour for small odd N destroys a fast convergence. The 
convergence of the gap sequences is slowed down by much stronger logarithmic corrections 
to the power-law behaviour than present in the ground state energy sequences. Thus the 
gained accuracy for higher spin quantum numbers is rather limited and larger rings together 
with methods like DMRG have to be used, see e.g. |3~1"|,|3"T|,|3"%| . 
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TABLE I. Properties of ground and first excited state of AF Heisenberg rings for various N 
and s: ground state energy Eq, gap AE, degeneracy deg, total spin S and shift quantum number 
k. f - O. Waldmann, private communication, ff - projection method |$(J|. Values for higher ./V 
are available from the author. 
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FIGURES 
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FIG. 1. L.h.s: Ground state energies (symbols) for antiferromagnetic coupled Heisenberg rings 
of s = 1/2 compared to the large N limit of Bethe and Hulthen (solid line). R.h.s: Relative 
deviation of the Levin u-estimate from the limit of Bethe and Hulthen. Plus symbols are used for 
even N, crosses for odd N. The circled symbols show how much the ground state energies itself 
deviate from the limit. 
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FIG. 2. L.h.s: Ground state energies (symbols) for antiferromagnetic coupled Heisenberg rings 
compared to the Levin ^-estimate (hatched area, sometimes shrunk to a line). Plus symbols are 
used for even N, crosses for odd N. 
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FIG. 3. Excitation energy of the first excited state for antiferromagnetic coupled Heisenberg 
rings and the Levin n-estimate (hatched area). Plus symbols are used for even N, crosses for odd 
N. To generate these figures also previously obtained exact diagonalization results of other authors 
have been used [ 26-29 , 31 1 . 
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